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Abstract 

A  multi  block  zonal  algorithm  which  solves  the  thin- 
layer  Navier-Stokes  and  the  Euler  equations  is  used  to  nu¬ 
merically  simulate  the  formation  and  roll-up  of  the  tip  vor¬ 
tex  in  both  subsonic  and  transonic  flows.  Four  test  cases 
which  used  small  and  large  aspect  ratio  wings  have  been 
considered  to  examine  the  influence  of  the  tip-cap  shape, 
the  tip-planform  and  the  free-stream  Mach  number.  It  ap¬ 
pears  that  both  the  tip-planform  and  the  tip-cap  shape  have 
some  influence  on  the  formation  of  the  tip  vortex,  but  its 
subsequent  roll-up  seems  to  be  more  influenced  by  the  tip- 
planform  shape.  In  general,  a  good  definition  of  the  forma¬ 
tion  and  the  roll-up  of  the  tip  vortex  has  been  observed  for 
all  the  cases  considered  here.  Comparisons  of  the  numerical 
results  with  the  limited,  available  experimental  data  show 
good  agreement  with  both  the  surface  pressures  and  the 
tip-vortex  strength. 

■Nomenclature 

B  =  semi  span  of  the  wing 

C  =  root  chord  of  the  wing 

Cp  -  drag  coefficient 

CL  =  lift  coefficient 


CP 

E,F,G 


quarter-chord,  pitching-moment  coefficent 
pressure  coefficient 
flux  vectors 

free-stream  Mach  number 
flow-field  vector 
velocity  vector 
Reynolds  number 
viscous  flux  vector 
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Uoo  =  free-stream  velocity 

X,Y,Z  =  Cartesian  coordinates 

a  =  angle  of  attack,  degrees 

Tv  =  dimensionless  strength  of  tip  vortex 

£ ,  r),  f ,  r  —  generalized  curvilinear  coordinates 

Introduction 

The  phenomenon  of  formation  and  subsequent  roll-up  of 
tip  vortices  behind  finite  aspect-ratio  wings  has  long  at¬ 
tracted  a  great  amount  of  attention  because  of  their  po¬ 
tential  hazard  to  aircraft  that  encounter  them  in  flight.  To 
date  there  are  many  theoretical,  numerical,  and  experimen¬ 
tal  studies  in  the  literature  which  are  devoted  entirely  to 
the  understanding  of  such  flows  encompassing  large  air¬ 
craft  wakes  (see  for  example  Refs.  1-2).  In  spite  of  this, 
the  present  understanding  of  such  flows  remains  essentially 
qualitative.  In  all  these  studies,  the  detailed  knowledge  of 
the  mechanics  of  the  formation  and  subsequent  roll-up  of 
tip  vortices  have  been  precluded.  In  addition  to  the  fixed- 
wing  aircraft,  the  operating  characteristics  of  helicopters 
are  strongly  influenced  by  the  vortex  wakes  of  the  rotat¬ 
ing  blades.  The  dominant  feature  of  such  wakes  is  the 
helical  vortices  which  have  their  origin  at  the  tips  of  the 
rotor  blades.  The  interaction  of  these  blades  with  the  vor¬ 
tex  wake  changes  their  aerodynamic  loads  and  thus  affects 
their  operating  performance,  vibration,  and  acoustic  char¬ 
acteristics.  Although  a  detailed  knowledge  of  the  formation 
and  the  initial  roll-up  of  such  concentrated  vortices  is  not 
used  for  typical  current  wake  analyses,  such  a  knowledge  is 
essential  to  modify  the  structure  of  these  vortices  and  to 
minimize  their  influence.  Thus,  it  is  important  to  study  in 
detail  the  flow  field  in  the  vicinity  of  the  tip  of  a  finite  wing 
and  of  a  helicopter  rotor  blade.  In  addition,  the  recent  in¬ 
terest  in  the  generation  of  high  lift  coefficients  for  V/STOL 
operations  has  spurred  further  increased  interest  in  under¬ 
standing  the  physics  of  the  formation  of  tip  vortices  and 
their  subsequent  downstream  trajectories. 

The  typical  flow  field  in  the  tip  region  of  a  finite  aspect 
ratio  wing  is  strongly  three-dimensional  and  often  is  sepa¬ 
rated.  Vortex  filaments  shed  in  this  region,  from  the  viscous 
layer  near  the  surface,  interact  with  each  other  to  form  a 


\t!\  1 


K 


discrete  tip  vortex  whose  strength  and  size  are  dictated  by 
the  aerodynamic  and  geometrical  constraints  of  the  wing 
and  of  the  wing  tip-geometry  respectively.  In  most  predic¬ 
tion  scheme-,  the  true  physical  situation  is  related  to  an 
analytical  model  through  the  concept  of  bound  vorticity. 
In  reality,  the  bound  vorticity  is  that  which  is  continually 
generated  a-  a  result  of  the  viscous  nature  of  the  bound¬ 
ary  layer  near  the  wing  surface  and  which  appears  to  be 
attached  to  the  wing  tip  in  the  Eulrrian  reference  frame 
(Kef. 

Therefore  it  appears  that  the  formation  of  tin  vortex 
originating  from  round  edge-  is  a  iscous  phenomenon. 
However.  tt>e  recent  work-  ot  Ric/t  et  al.  (lie:.  1  and 
of  Mitr  heltree  et  al  (Ref.  5j  have  used  Euler  equation, 
to  numerically  simulate  the  vortex  formation  and  its  sun- 
sequent  roll-up.  The  reason  for  the  success  of  such  a  pro- 
cedute  i-  that  the  artificial  viscosity,  used  in  the  nunerical 
-i  heme  to  <•<■!. trol  the  nonlinear  stability,  has  piayed  the 
role  o!  the  na  ural  vi.tosity.  More  lecently.  in  a  pioneering, 
hut  iimited  work.  Mansour  (Ref.  (i)  used  the  thin-layer 
Navier-stokes  equations  (lief  7  .  in  conjunction  with  an 
elliptic  grid  generation  scheme,  to  .  itnulate  the  flow  held 
anti  the  tin-vortex  of  a  low -aspect-ratio  wing.  In  spite  of 
a  (lever  data  management  strategy,  the  central  p'ot es-itig 
unit  (t  PI  )  time  lot  each  -tt  ady-state  calculation  exceeded 
2‘  hours  ot  NASA  Ames  CRAY  1-S  machine,  veil,  only 
7'..00(i  grid  points  in  the  entire  Ho  >  field,  for  a  fivc-order- 
ot-magnit uth  dtot>  it  > !’«*  residuals.  In  the  pre-em  work 
an  attentate  but  n, ore  eflicient  met  he  d  of  multi  block  zo  ia’ 
approach  (Rtfs.  >•.!•)  is  used  to  simulate  the  Tow  field- 
of  isolated  wing  sections  by  solving  the  Euhr  thin-layer 
Navier-Stoke?  etptations  with  a  view  to  simulating  the  tip- 
vortex  formation  in  subsonic  and  transonic  flows. 


Governing  Equations  and  Numerical  Scheme 


The  governing  partial  differential  equations  are  the  thin- 
layer  Navier-Stokes  and  the  Euler  equations  (Refs.  7.10). 
The  equations  are  generally  nondimensionalized  by  frec- 
stream  quantities  and  are  transformed  to  the  arbitrary 
curvilinear  spare  (J.  q.  y.  r)  while  retaining  strong  conser¬ 
vation  law  form  to  capture  shock  waves.  The  transformed 
equations  written  in  generalized  rurvilineai  coordinates  art 
given  by  (Refs.  7.10) 
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At  the  Rey  nolds  numbers  of  it. 'ores  for  the  preset!  cal¬ 
culation?.  tl  e  flow  around  ti  e  wing  ca  t  be  considered  to  be 
mostly  itivi-cid  except  in  a  small,  bit'  crucial,  region  near 
the  body  si  rface  when  the  viscous  effects  are  important. 
Accordingly,  the  flow  held  is  solved  using  a  multi  block 
zonal  approach  developed  at  NASA  An.es  Research  Ou¬ 
ter  (Refs,  s.'.i)  wherein  the  invisrid  outer  blocks  are  solved 
using  'he  Euh-r  equations  and  the  inner  viscous  blocks  ate 
solved  using  the  thin-layer  Navier-Stokes  equations.  The 
results  of  w  ing  flow  topology  and  of  the  tip-vortex  are  pre¬ 
sented  for  four  different  wing  shapes,  with  three  of  them 
being  in  tran-onic  flow  and  one  in  subsonic  flow.  The  in¬ 
fluences  of  the  tip-planform  and  of  the  tip-cap  effect  on 
the  tip-vortex  formation  are  also  discussed.  In  addition, 
the  results  from  the  multi  block  zonal  approach  are  com¬ 
pared  with  the  results  from  a  one-block  solution  of  thin- 
layer  Navier-stokes  equations.  The  numerical  results  are 
compared  with  the  limited  available  experimental  data. 


and  (  0  for  Euler  equations,  and  <  -  1  thii-laver 

Navier-Stoke?  equations.  The  viscous  flt.x  vector  n  the 
thin-layer  approximation,  is  given  in  Re's.  /  at  d  ’0.  The 
primitive  variables  of  Eq.  (1).  viz.,  p, pu.pr.pt/'  and  e  are 
normalized  by  the  free-stream  reie  ence  quantities.  The 
rein t ion?  for  the  contravariat  t  velocities  ai  d  IV.  the 

Jacobian  of  transformation  J,  and  tie  met  ic.  are  lescrihed 
in  Refs.  T  and  10. 

The  velocity  components  u.v.w  and  the  pressure,  p.  are 
related  to  the  total  energy  per  unit  volume,  c.  through  the 
equation  of  state  for  a  perfect  gas  by 

P  (‘  D(r  £(u-  •  r-  •  (('*'))  (it) 

The  numerical  code  which  solves  these  equations,  called 
AROl)  (Refs.  10).  was  adapted  by  Holst  et  al  (Ref  *) 
and  Flores  (Ret.  **J  to  develop  the  transonic  Navier-Stoke- 
(TNS)  multi  block  zonal  algorithm  Most  of  the  impor¬ 
tant  features  of  the  AR('dl)  rode  are  eontdi.cd  in  th  *  TN- 
rode  with  one  important  difference,  viz.,  w  hen  ilea  ling  w  it  n 
multiple  blocks.  The  development  of  this  4-hlcxk  version  of 
the  ’TNS  code  and  tur  improvements  to  the  numerical  al¬ 
gorithm  to  handle  multiple  zones  or  block?  are  described 


in  detail  in  the  papers  of  Holst  et  al.  (Ref.  8)  and  of 
Flores  (Ref.  9).  In  addition  to  the  example  solutions  for 
isolated  wings  discussed  in  these  two  papers.  Kaynak  et  al. 
(Ref.  11)  report  further  refinements  to  this  rode  to  analy  ze 
large-scale  separation  on  wings.  The  present  version  of  this 
4-block  scheme  is  mainly  used  for  computing  isolated  wing 
solutions,  with  or  without  wind-tunnel  walls.  However,  a 
16-block  version  of  this  code  for  simulating  the  complete 
flow  field  of  a  fighter  aircraft  is  currently  under  develop¬ 
ment. 

There  are  five  important  features  of  the  TNS  code.  (I) 
The  thin-layer  Navier-Stokes  equations  and  the  Euler  equa¬ 
tions  are  solved  in  strong  conservation-law  form  to  cap¬ 
ture  the  shock  waves.  (2)  The  convergence  procedure  is 
significantly  accelerated  over  the  standard  ARC3D  code 
(Ref.  10).  which  solves  thin-layer  Navier-Stokes  equations 
over  the  complete  flow  domain,  because  the  present  scheme 
solves  the  Euler  equations  over  a  significant  part  of  the  flow 
domain.  (3)  There  are  two  numerical  options  for  solving 
the  equations.  One  is  based  on  the  ADI  algorithm  which 
solves  block-tridiagonal  matrices  along  each  coordinate  di¬ 
rection  (Ref.  12).  The  other  is  based  on  the  diagonalized 
algorithm  (Ref.  13)  which  solves  a  set  of  five  scalar  penta- 
diagonal  matrices  along  each  coordinate  direction.  (4)  The 
diagonal  algorithm  has  been  implemented  with  two  options 
of  variable  time-step  philosophies  to  accelerate  the  conver¬ 
gence  rate  of  the  numerical  scheme,  viz.,  one  that  scales 
the  marching  time-step  with  the  transformation  Jacobian 
(Ref.  14)  and  the  other  that  scales  the  time-step  using  a 
combination  of  the  Jacobian  and  local  solution  variation 
(Ref.  9.10).  (">)  The  code  is  vectorized  for  the  C'RAY-XMP 
machine. 

Both  numerical  schemes  (  TNS  and  AIH'3!)  )  use  the 
standard  see  ond-order-accurate  central  differencing  to  con¬ 
struct  the  appropriate  spatial  differencing  scheme.  The 
diagonal  scheme,  which  is  used  in  the  present  study  to 
calc  ulate  steady-state  solutions,  uses  fourth-order-accurate 
smoothing  operators  on  the  implicit  and  explicit  'ides  of 
the  numerical  algorithm  tor  controlling  nonlinear  stability 
of  the  numerical  scheme  In  the  present  c  ale  ulat  ions  a  tur¬ 
bulent  boundary  layer  i-  assumed  for  the  entire  wing,  and 
the'  Baldwin  and  Fornax  algebraic  turbulence  model  (Ref. 
III.  which  i-  currenth  available  in  the  TNs  code,  (s  used 
to  calculate  the  turbulent  eddv  vi-co-m 
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is  shown.  The  coarse  (inviscid)  outer  grid,  marked  as  block 
(1).  is  shown  in  white,  the  finer  inviscid  grid,  marked  as 
block  (2).  is  shown  in  red.  The  two-blocks  adjacent  to  the 
wing  and  on  either  side  of  it.  marked  as  blocks  (3)  and  (4| 
are  shown  in  yellow  .  and  are  the  viscous  blocks  which  have 
clustering  in  the  leading-edge  and  t railing-edge  regions  as 
well  as  in  the  wing  tip  region  which  is  used  to  resolve  the 
tip  vortex.  The  inviscid  grids  also  have  clustering  as  shown. 

In  the  present  case,  each  of  the  four  zones  ty  pically  has 
over  40.000  grid  points.  Although  the  grid  geometry  for 
each  wing  discussed  here  depended  on  the  aspect  ratio  of 
the  individual  wing,  the  spanwise  grid  spacing  at  the  wing 
tip  and  the  spacing  in  the  normal  direction  to  the  wing 
surface,  which  were  used  to  resolve  the  boundary  layer,  were 
kept  same  for  all  the  wings,  viz..  0.015  of  chord  and  1  x 
10  5  of  chord,  respectively.  The  viscous  blocks  had  25  grid 
points  in  the  normal  direction  for  all  the  wings. 

The  data  managemant  strategy  is  discussed  by  Holst  et 
al.  (Ref.  8)  and  also  by  Kaynak  et  al.  (Ref.  11).  In  brief, 
the  base  grid,  usually  generated  outside  the  TNS  code,  is 
“read  in"  first  and  this  grid  is  divided  into  the  proper  zones 
by  the  ‘"Zoner"  code.  Once  this  is  done,  the  flow  solver  is 
initiated.  The  iteration  procedure  starts  in  the  outer  in¬ 
viscid  zone  or  block,  and  proceeds  into  the  inner  viscous 
blocks.  The  information  necessary  to  update  the  boundary 
conditions  at  the  zonal  interface  is  found  from  the  neigh¬ 
boring  zones  through  a  series  of  one-dimensional  linear  in¬ 
terpolations.  Such  a  scheme  lends  itself  to  a  conservative 
treatment  of  the  boundaries  and  thus  captures  distortion- 
free-movement  of  discontinuities  across  the  boundaries. 

In  the  solution  procedure,  only  the  information  neces¬ 
sary  to  solve  each  zone  resides  at  any  one  time  in  the  main 
memory  of  the  CRAY-X.MP  computer.  The  information  of 
zones  which  are  not  being  computed  is  stored  temporarily 
on  the  solid  state  device  (SSD).  The  use  of  SSD  allows  a 
great  deal  of  flexibility  and  thus  frees  the  main  memory  of 
the  spare  normally  taken  by  the  data  stored  on  the  SSD 
and  enables  the  use  of  maximum  number  of  grid  points  in 
the  flow  field.  Also,  the  use  of  SSD  reduces  the  I  O  wait 
time  significantly  (Ref.  9). 

Houndary  Conditions 

Since  the  TNS  rode  is  a  multi  zonal  algorithm,  there  are 
two  t v  pe-  of  boundaries  where  conditions  have  to  be  spec¬ 
ified  viz..  (II  the  physical  boundaries  such  as  inflow,  out¬ 
flow.  and  the  surface  conditions:  and  (2)  the  zonal  bound¬ 
aries  where  all  the  flow  quantities  have  to  be  made  contin¬ 
uous  V 1 1  the  boundarv  and  interface  or  zonal  conditions 
used  are  applied  explicitly 

\t  tin’  far-held  boundaries,  which  are  6-10  chords  away 
from  the  vvitig.  Iree-strearn  values  are  specified  for  all  the 
flow  quantities  \t  (tie  outflow  boundaries,  zeroth-order 
extrapolation  o  used  from  the  grid  interior  \t  the  sym- 
rnetrv  plane,  a  zeroth-order  extrapolation  is  used  for  the 
densitv  and  a  first-order  extrapolation  is  used  for  the  x- 
compoiicnt  and  the  /-component  of  velocities  while  setting 


the  spanwise  velocity  component  to  zero  to  force  symme¬ 
try.  A  first-order  extrapolation  is  also  used  for  the  pressure 
and  the  energy  is  calculated  from  the  equation  of  state.  At 
the  surface  of  the  wing,  a  no-slip  condition  is  used  for  the 
three  velocity  components  and  the  pressure  is  calculated  by 
solving  the  normal  momentum  equation  at  the  surface.  In 
the  spirit  of  the  thin-layer  approximation,  the  normal  mo¬ 
mentum  equation  would  translate  into  dp  dn  0.  where  n 
is  the  local  normal  to  the  surface.  Density  is  determined 
by  assuming  adiabatic  wall  conditions.  Having  known  the 
pressure,  p.  and  density,  p.  at  the  wall,  the  total  energy,  e, 
is  determined  from  Kq.  (3).  The  details  of  the  interface 
boundary  conditions  at  the  zone  interfaces  are  described  in 
Refs.  8  and  1 1 . 

Results  and  Discussion 

In  this  section,  numerical  results  for  flow  fields  and  tip- 
vortex  formation  are  presented  for  four  different  wing  sec¬ 
tions.  These  results  are  compared  with  the  available  exper¬ 
imental  data.  The  different  flow  conditions  and  the  wing 
geometries  considered  consist  of:  (1)  a  rectangular  wing 
with  an  aspect  ratio  of  2.5  without  twist  or  taper,  in  a  uni¬ 
form  flow  of  A =  0.16  and  at  11  degrees  angle  of  attack; 
(2)  a  wing  with  a  small  aspect  ratio  of  0.83  with  twist  and 
taper  in  a  uniform  flow  of  Mx  =  0.9  and  at  5  degrees  angle 
of  attack;  and  (3)  two  wings  with  exotic  tip  shapes,  both 
having  an  aspect  ratio  of  5.0,  in  a  uniform  flow  of  = 
0.85  and  at  5  degrees  angle  of  attack.  The  planforms  and 
the  fine  surface  grid  (at  every  grid  point)  for  these  different 
wing  geometries  are  shown  in  Fig.  2.  The  boundary  layer 
is  assumed  to  be  turbulent  for  the  entire  wing  sections. 

Rectangular  Wing  in  Subsonic  Flow 

This  wing  consists  of  a  rectangular  wing  with  a  square 
tip  and  is  made  up  of  NACA-0015  airfoil  sections  without 
twist  or  taper.  Its  planform  and  surface  grid  are  shown  in 
Fig.  2a.  For  this  wing  in  a  uniform  free-stream  of  .V/oc  = 
0.16  and  at  1 1  degrees  angle  of  attack,  the  steady-state  flow 
field  was  calculated  using  the  T.NS  code.  The  flow  Reynolds 
number,  based  on  the  free-stream  velocity  and  the  chord  of 
the  wing,  was  2  million.  Typical  computational  time,  for  a 
four-order  drop  in  residuals,  was  of  the  order  of  3  hours  on 
the  CRAY-XMP  machine. 

The  steady -state  results  for  this  w  ing  are  show  n  in  Figs. 
3-6.  Figure  3  shows  the  computed  surface  pressure  dis¬ 
tributions  at  several  spanwise  stations  compared  with  the 
experimental  data  of  Spivey  and  Moorhouse  (Refs.  16.17). 
The  inboard  spanwise  stations  show  good  agreement  with 
the  experiments  as  seen  in  Figs.  3a-3c.  However,  the  suc¬ 
tion  peak  associated  with  the  tip  vortex  that  was  observed 
at  the  97%  spanwise  station  near  the  trailing  edge  region, 
as  shown  in  Fig.  3d.  has  not  been  well  predicted  The 
probable  reason  for  this  ran  be  explained  by  examining  the 
tip-cap  profile  of  the  present  wing  goemetry.  The  top  half 
of  Fig.  4  shows  the  details  of  the  tip-cap  at  the  X  50% 
station.  This  tip-cap.  although  it  resembles  a  circular  cap 
at  this  X-station.  it  has  a  very  different  shape  at  other 
stream  wise  stations.  A  detailed  examination  of  the  tip- 


cap  profile  at  other  streamwise  stations,  shown  for  example 
for  X  =  75%  and  for  X  =  90%  in  Fig.  4,  indicates  that  its 
shape  changes  in  the  leading-  and  trailing-edge  regions  to 
that  of  a  sharp  corner.  This  sharp  corner  at  the  leading- 
edge  apparently  may  be  responsible  for  producing  a  larger 
separation  than  a  square-tip  wing  would  produce  and  thus 
nullifying  the  suction  peak  associated  with  the  tip  vortex. 
The  calculated  pressure  field  agrees  with  the  experiment  at 
this  97%  spanwise  station  on  the  lower  surface  and  on  most 
of  the  upper  surface,  except  for  the  important  suction  peak 
near  the  trailing-edge  region. 

To  examine  further  the  influence  of  the  tip-cap  shape 
on  the  flow  field  in  the  tip  region,  the  rectangular  wing 
with  square  tip  was  modified  by  adding  an  ellipsoidal  tip- 
cap  that  produces  a  smooth,  rounded  surface  along  the  tip. 
This  cap  has  a  minor  axis  equal  to  half  the  thickness  of 
the  wing.  In  order  to  represent  this  smooth  curved  tip.  the 
surface  grid  resolution  was  increased  to  twice  the  value  in 
the  spanwise  direction  to  that  used  earlier.  The  tip  of  the 
wing  with  this  tip-cap  is  rounded  off  in  all  the  directions 
to  avoid  sharp  corners.  The  finite  difference  grid  for  this 
wing  consists  of  a  spherically  warped  grid  of  0-0  topol¬ 
ogy  generated  by  the  hyperbolic  grid  solver  of  Steger  and 
Chaussee  (Ref.  18).  and  is  shown  in  Fig.  4  at  the  X  -  50% 
streamwise  station. 

To  generate  the  steady-state  solution  for  this  wing  at  the 
same  flow  conditions  of  Mx  -  0.16.  a  11  degrees,  and 
Re  =  2  million,  a  single-block,  thin-layer  Xavier- Stokes  flow 
solver  called  PEN3D  (Ref.  19)  is  used.  This  code  is  also  a 
derivative  of  the  ARC3D  code  (Ref.  10)  but  is  coded  in  the 
VECTORAL  language.  The  essential  details  of  this  code 
are  described  in  Pan  and  Pulliam's  paper  (Ref.  19).  All  the 
features  of  this  code  are  identical  to  the  TNS  code  but  it 
is  limited  to  applications  with  only  single-block  grids.  The 
present  single-  block  grid  has  the  dimensions  of  110  points 
in  the  periodic  direction.  44  in  the  spanwise  direction,  and 
66  in  the  normal  direction. 

Typical  results  of  the  steady-state  solution  for  this  wing 
are  shown  in  Figs.  3c-d  in  the  form  of  surfare-pressure  dis¬ 
tributions.  The  results  agree  well  with  the  TNS  calculations 
and  the  experimental  data  of  Spivey  and  Moorhouse  (Refs. 
16.17)  from  the  wing  root  to  the  90%  spanwise  station.  Al¬ 
though  the  PEN3D  calculations  show  a  suction  peak  in  the 
tip  region,  because  of  the  formation  of  the  vortex,  it  is  not 
as  large  as  that  observed  in  the  experiment 

This  rounded  tip-cap  did  not  produce  any  separation  in 
the  tip  region  on  the  upper  surfaceof  the  wing,  as  seen  from 
the  surfare  oil  picture  of  Fig  5a  (The  surface  oil  flow  pic¬ 
ture  is  generated  b\  releasing  fictitious  fluid  particles  at  one 
grid  point  above  the  surface  and  by  restricting  these  parti¬ 
cles  to  that  plane  1  In  c  out  rast .  t  he  t  ip-e  ap  shape  discussed 
earlier  produce-  large  separation  in  the  tip  region  for  this 
wing  as  shown  in  the  -urface  oil  How  pic  ture  of  Fig  6  The 
extent  ol  the  separation  region  on  this  wing  is  more  pro¬ 
nounced  than  i-  apparent  from  the  pressure  contour  plot  of 
Chigier  and  (  or-igl.a  (Ref  20)  indicating  that  the  tip-cap 
of  Fig  4a  top  \  lew  .  tir-i  thought  to  mimic  the  experimental 


square  tip.  is  really  different  and  is  neither  a  near-square 
nor  a  near-rounded  tip.  The  sharp  corner  at  the  leading 
edge  in  the  tip  seems  to  he  separating  the  flow  at  this  point 
as  seen  from  the  close-up  view  of  surface  oil  picture  of  Fig. 
7a. 

The  close-up  views  of  the  lift-off  of  the  tip  vortex  for  these 
two  tips  are  shown  in  Figs.  5b  and  7b.  The  tip  vortex  for 
the  rounded  tip  wing  of  Fig.  5a  is  lifting  off  from  the  tip- 
face  well  before  the  trailing  edge.  For  the  non-rounded  tip 
of  Fig  7b.  the  tip  vortex  is  lifting  off  from  the  upper  surface 
of  the  wing  at  about  the  trailing  edge  position.  Both  the 
tip  shapes  of  Figs.  5b  and  7b  show  the  three-dimensional 
trajectories  of  the  particles  around  the  tip.  The  particles 
released  on  the  lower  surface  of  the  wing  cross  over  the  tip 
region  into  the  low-pressure  region  of  the  upper  surface  of 
the  wing.  These  particles  mix  with  the  particles  released 
on  the  upper  surface  and  together  they  define  a  tip  vortex 
that  is  distinct  from  the  rest  of  the  sheet.  The  tip  vortex 
then  rolls  up.  slightly  inboard  of  the  wing  tip.  as  it  leaves 
the  wing  surface  (Fig.  5b).  Further  downstream  of  this,  the 
tip  vortex  rolls  inboard  and  stays  distinctly  above  the  shed 
vortex  sheet  as  seen  in  the  far-field  view  of  the  tip  vortex 
shown  in  Fig.  8.  Also,  shown  in  Fig.  *  are  the  vorticity 
contours  in  the  Y-Z  planes  through  the  tip  vortex  at  sev¬ 
eral  downstream  X-locations.  In  this  figure,  the  magnitude 
of  the  vortex  contour  increases  in  the  order  of  the  change 
in  color  according  to  the  blue-green-yellow-red.  Judging 
from  the  magnitudes  of  the  contours,  the  core  vorticity  is 
maximum  when  it  is  just  downstream  of  the  trailing  edge 
of  the  wing.  However,  at  the  downstream  grid  boundary, 
which  is  6  chords  away  from  the  wing,  the  vorticity  levels 
have  reduced  by  an  order-of-magnitude  caused  by  a  very 
coarse  grid  structure.  It  is  also  interesting  to  see  that  the 
vortex  structure  in  the  far-field  is  stretched  in  accordance 
with  the  local  grid  shape.  Close  to  the  wing,  the  structure 
is  somewhat  axisymmetric  and  as  the  grid  gets  coarser  in 
block  I.  the  vortex  shape  gets  stretched  accordingly 

The  lift.  drag,  and  quarter-chord  pitching-moment  coeffi- 
cents  calculated  from  the  TNS  code  for  this  wing  are  0.702. 
0.0748.  and  -0.00076  respectively.  This  compares  well  w  ith 
the  measured  values  of  0.763.  0.0868.  and  -0.00813  for  the 
lift,  drag,  and  pitching-moment  coefficients,  respectively, 
at  o  12  degrees.  The  lift  per  degree  of  angle  of  at¬ 
tack  .  ('i  o.  for  both  the  experiment  and  the  TNS  cal¬ 
culations  is  approximately  0.064  degree.  The  modified  tip 
with  PE.N.'if)  gave  a  value  of  0.07  degree.  From  this  the 
approximate  strength  of  tip  vortex,  using  the  definition  of 
lift  equivalent  vortex  strength,  is  0.35.  (Note  the  vortex 
strength  quoted  here  is  a  dimensionless  quantity,  nondi- 
mensionli/ed  by  the  free-stream  velocity  and  by  the  chord 
of  the  wing  |  Another  method  of  evaluating  the  tip  vortex 
strength  i-  In  doing  a  line  integral  of  the  velocity  vector 
over  a  <  lo-ed  path  enclosing  the  vortex,  viz.. 
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The  result  of  such  a  line-integral  at  several  X-stations 
behind  the  wing  over  a  path  big  enough  to  enclose  the  vor¬ 
tex  gave  a  range  of  values  from  0.33-0.38.  depending  on  the 
size  of  the  integration  path  around  the  vortex.  The  larger 
value  corresponds  to  the  case  of  line-integral  pal  h  exending 
all  the  way  to  the  wing  root  in  the  ^'-direction  and  extend¬ 
ing  to  the  limits  in  Z-direction.  This  path  will  include  the 
contribution  from  the  vortex  sheet  in  the  wake  also.  The 
t if»  vortex  strength  from  these  two  methods  is  estimated  to 
be  0.35. 

Wing  C  in  Transonic  Flow 

W  INC;  C  belongs  to  the  class  of  advanced  technology 
wing  sections  and  has  been  extensively  studied  both  com¬ 
putationally  (Refs.  6.11)  and  in  wind-tunnel  experiments 
(Refs.  21.22).  Its  planform  and  the  surface  grid  are  shown 
in  Fig.  2b.  It  is  a  low-aspect-ratio  wing  (0.83  based  on  the 
root  chord)  which  is  made  of  supercritical  wing  sections, 
with  a  twist  angle  of  8.17  degrees,  a  taper  ratio  of  0.3  and 
a  leading-  edge  sweep  of  45  degrees.  Recently  Kaynak  et 
al.  (Ref.  11)  have  presented  extensive  computational  re¬ 
sults  for  this  wing  for  a  range  of  Mach  numbers.  These 
results  were  also  computed  using  the  TNS  rode.  Whereas 
this  study  concentrated  on  analysing  the  flow  topology,  the 
present  sludv  specializes  in  the  flow  in  the  vicinity  of  the 
wing  tip  and  in  particular  concentrates  on  the  tip  vortex 
phenomenon.  As  noted  earlier.  Mansour  (Ref.  6).  in  a  pi¬ 
oneering  but  limited  study,  attempted  to  simulate  the  tip 
vortex  formation  for  this  wing  in  a  free-stream  Mach  num¬ 
ber  of  0.82  and  at  5  degrees  angle  of  attack. 

The  flow  conditions  for  the  present  computations  are  A/v 

0.9.  o  5.  and  Re  6.8  million  based  on  the  mean 
aerodynamic  chord.  As  before,  the  baseline  coarse  grid  was 
generated  by  the  parabolic  grid  solver  of  Edwards  (Ref.  15) 
with  sufficient  grid  resolution  in  the  wing  tip  region  so  as 
to  resolve  the  tip  vortex.  The  interior  grids  are  generated 
within  the  TNS  code  by  the  Z oner  program.  The  CPI' 
time  for  a  fully  converged  solution  was  about  3  hours  for 
this  case  which  is  highly  transonic  and  has  a  shock-induced 
flow  reparation  on  the  upper  surface  of  the  wing. 

The  calculated  surface  pressure  distributions  for  several 
s panwise  stations  are  shown  in  Fig.  9  and  are  compared 
wit  b  t  he  experimental  data  of  Keener  (Ref.  21 ).  The  results 
show  good  agreement  with  the  experiments  over  the  parts 
of  the  wing  which  do  riot  have  massive  shock-induced  sep¬ 
aration  However,  the  --hock  wave  position  and  post-shock 
pressure  levels  are  not  well  predicted  at  the  90'<  span  sta¬ 
tion.  Similar  results  were  obtained  by  Kay  nak  et  al.  (Ref. 

1 1 ).  w  bo  discuss  t  lie  possible  reasons  for  this,  including  in¬ 
adequate  grid  resolution,  artificial  dissipation,  and  turbu¬ 
lence'  model.  The  other  features  of  the-  flow  are  the  presence 
of  a  strong  shock  wave  in  the  leading  -edge  region  extend¬ 
ing  all  t hr  way  into  the  inboard  of  the  tip  region.  Here  it 
coalesces  with  the  secondary  shock .  This  combined  shock 
produces  large  boundary -layer  separation.  Figures  lOa-b 
show  this  shock  at  'Ml1';  spamvise  station,  and  llu  extent  of 
separation  is  shown  in  Fig.  l()c.  This  figure  is  generated 
by  releasing  the  fluid  particles  at  one  grid  point  above  the 


surface  and  by  not  restricting  them  to  any  surface.  The 
curling  of  the  fluid  particles  in  the  separated  region  indi¬ 
cates  the  presence  of  vortical  flow.  Figures  11-12  show  for 
the  upper  surface  of  the  wing  the  surface  oil  flow  pattern 
and  the  velocity  vector  plot.  (Just  as  in  the  surface  oil 
picture,  the  velocity  vector  plot  is  generated  by  releasing 
particles  at  one  grid  point  above  the  surface).  Both  these 
figures  clearly  show  the  extent  of  separation. 

Tip-vortex  formation  was  visualized  by  releasing  several 
particles  at  different  locations  along  the  span  and  various 
heights  from  the  wing  surface.  Figures  1 1-12  show  the  close- 
up  and  far-field  view  of  the  tip  vortex.  Because  of  the  large 
cross  flow  caused  by  the  high  sweep  angle,  the  boundary 
layer  in  the  tip  region  becomes  thicker.  This  thick  bound¬ 
ary  layer  easily  separates  in  the  presence  of  the  shock  wave. 
Fart  of  this  separated  How  lifts  up  from  the  wing  surface 
and  merges  into  the  flow  in  the  tip  region.  These  fluid 
particles  get  entangled  with  the  particles  coming  over  the 
tip  from  the  high-pressure  region  of  the  wing's  lower  sur¬ 
face  and  in  combination  these  particles  separate  into  the  tip 
vortex,  and  are  distinct  from  the  rest  of  the  vortex  sheet. 
This  action  is  apparent  from  the  picture  of  Fig.  Ill  which 
shows  the  close-up  view  of  the  formation  process  of  the  tip 
vortex.  The  lift-up  of  the  tip  vortex  from  the  surface  is 
apparent  in  this  view.  Figure  1-1  shows  the  far-field  view  of 
the  same  vortex  .  This  view  shows  how  the  t ip  vortex  rolls 
inboard  in  the  downstream  wake.  The  vorticity  contours  in 
the  Y-Z  planes  at  several  X-locations  are  seen  in  fig.  1-4. 
As  noted  before,  the  magnitude  of  vorticity  decreases  in 
the  far-field  coarse  grid.  But  the  line-integral  of  the  veloc¬ 
ity  vector  around  a  dosed  path  surrounding  the  tip  vortex 
would  give  the  same  value  for  the  strength  of  vortex  as  long 
as  the  integration  path  was  large.  For  this  wing,  the  cal¬ 
culated  values  of  the  aerodynamic  coefficients  were 
0.301.  (  0.0317.  and  C\,  -0  110.  The  strength  of  the 

vortex  determined  from  the  line-integral  method  yielded  a 
value  between  0.17  and  0.185.  depending  on  the  size  of  the 
integration  path  around  the  tip  vortex.  The  smaller  value 
includes  the  vortex  sheet  in  the  wake  and  therefore  corre¬ 
sponds  to  the  larger  path  of  t lie  integral  (usually  extending 
to  the  sy  linnet r>  plane  in  the  Y-directioti  and  limit'  of  Z  in 
the  Z-dirert  ion.  This  implies  that  the  tip  vortex  strength  is 
equal  to  0.1  No  (normalized  by  the  free-stream  velocity  and 
the  root  chord  of  the  wing). 

ONER  A  Wing  in  Transonh  Flow 

This  wing  has  an  aspect-ratio  of  3  and  represents  a  typ¬ 
ical  helicopter  rotor  blade.  The  tip  of  the  wing  is  more 
exotic  than  are  the  tips  of  the  two  previous  wings.  The 
wing  is  made  up  of  three  special  airfoil  sections  ranging  in 
thickness  from  12  to  0  ‘7.  the  (flicker  section  being  at  the 
root  arid  the  thinnest  one  at  the  tip  The  vvirig  Ini'  a  con¬ 
stant  chord  up  to  about  70’7  of  the  spanwise  station  and 
then  dei  reas.  -  m  t fu  ratio  of  the  thickness  ratio.  Consid¬ 
ering  that  wi  i  an  have  only  23  points  ori  the  wing  in  the 
spanwi'e  direition  llmiited  h\  the  memory  of  the  CRAY- 
XMIM  in  the  fine  filler  and  the  visions  grid',  the  same 
number  of  points  a-  we  had  for  the  WIMi  C.  the  spanwise 


grid  distribution  will  be  very  sparse  for  this  wing.  Figure 
2c  shows  the  planform  and  the  surface  grid  for  this  wing. 
For  this  wing  in  a  uniform  free-stream  of  .W-*  -  0.85.  a  = 
5  degrees  and  Re  8.5  million  (based  on  the  root  chord), 
computations  were  carried  out  to  generate  a  steady-state 
solution. 

Figure  15  shows  the  computed  surface-pressure  distibu- 
tions  for  several  spanwise  stations  along  the  wing.  Exami¬ 
nation  of  the  pressure  distributions  indicate  the  presence  of 
two  shock  waves,  a  weak  leading-edge  shock,  and  a  strong 
shock  downstream  of  this.  The  Mach  number  contour  plots 
of  Fig.  1C  confirm  this.  (The  Mach  number  contours  were 
drawn  at  a  plane  outside  the  viscous  region.)  The  main 
shock  extends  all  the  way  to  the  root  of  the  wing  which 
is  a  symmetry  plane.  This  shock  produces  large  scale  sep¬ 
aration  on  the  wing.  The  surface  oil  How  pattern  of  Fig. 
17  shoves  the  separation  and  reattachment  lines  along  with 
other  flow  features  clearly.  The  shock  becomes  weak  to¬ 
wards  the  wing  tip  and  the  flow  gets  complicated  between 
tlie  75‘7  and  90‘7  spanwise  stations.  An  outward  spiralling 
vortex  emanates  from  the  separated  flow  in  this  region. 

Figures  18-19  shows  two  views  of  the  tip  vortex  for  this 
wing.  The  initial  formation  of  the  vortex  and  lifting  off 
from  the  wing  surface  is  shown  in  Fig.  18.  It  appears  that 
the  swept-back  tip  has  similar  influence  as  the  sweep  effect 
(see  Wing  O  results  above),  viz.,  the  cross  flow  and  tip¬ 
loading  influence.  This  cross  flow  enables  the  fluid  particles 
released  in  the  vicinity  of  the  tip  to  first  braid  and  then  to 
roll-up  and  lift  off  from  the  surface.  Also,  the  fluid  particles 
released  on  the  high-pressure  side  of  the  wing  (lower  sur¬ 
face)  at  the  up.  cross  over  to  the  low-pressure  side  (upper 
side)  by  wrapping  around  the  wing  tip  in  the  classical  three- 
dimensional  tip-relief  fashion  and  braid  in  to  the  swirling 
tip  vortex.  Figure  19  shows  the  far-field  view  of  this  tip 
vortex.  The  lift  off  of  the  tip  vortex  from  the  wing  sur¬ 
face  and  then  t he  roll  inboard  in  the  downstream  wake  are 
clearly  ev  ident  from  this  photograph.  Also,  shown  in  this 
are  vorticity  contours  drawn  in  the  Y-Z  planes  through  the 
tip  vortex  at  several  X-lorations  behind  the  wing.  As  ob¬ 
served  before,  the  vorticity  levels  of  the  contours  decrease 
in  the  downstream  roarse-grid  region  of  Block  1.  and  the 
shapes  of  these  contours  get  stretched  in  the  pattern  of  the 
grid  geometry 

The  lift.  drag,  and  the  pitching-moment  coefficients  for 
this  wing  art1  f 0.20.  (’/>  0.061.  and  f’.\(  0.002. 

Note  that  the  pitching-moment  is  positive  (destabilizing) 
for  this  wing.  Based  on  the  definition  of  lift-equivalent  vor¬ 
tex  strengt  h.  the  non  dimensional  tip  vortex  strength  is  0.13. 
Also,  the  value  of  I\-  determined  from  the  line-integral  of 
the  velocity  vector  around  a  closed  path  surrounding  the 
vortex  and  at  several  X-lorations  in  the  wake  lies  in  the 
range  between  0.12  -  0.11  depending  on  the  path  of  the 
integral.  The  larger  value  corresponds  to  the  the  smaller 
integral  path  surrounding  the  tip  vortex  The  smaller  value 
includes  the  effect  of  wake  vortex  sheet.  From  this  calcu¬ 
lation  it  can  fie  rone  hided  that  the  tip  vortex  strength  is 
0  1  1 


Effect  of  Tip  Planform 

To  study  the  influence  of  the  tip  geometry  on  the  tip- 
vortex  formation,  the  ONERA  wing  of  Fig.  2c  was  modified 
to  reshape  the  planform  in  the  tip  region,  keeping  the  rest  of 
the  geometry  the  same.  The  planform  and  the  surface  grid 
for  this  modified  wing  is  shown  in  Fig.  2d.  This  modified 
wing  geometry  has  an  unswept  leading  edge  and  looks  more 
like  a  rectangular  wing  except  the  tip  region  retains  the 
original  taper.  The  flow  conditions  are  kept  identical  to 
the  case  of  ONERA  wing  of  section  (c)  above.  Steady-state 
solution  was  generated  as  before  using  the  TNS  code  for  this 
wing  section.  From  this,  the  lift,  drag,  and  the  pitching- 
moment  coefficients  for  this  wing  are  C;.  -  0.332,  Cp  — 
0.0796,  and  C*/  -  -0.018.  The  obvious  gross  changes  in 
flow  features  are  the  increase  in  the  lift  and  drag  values,  and 
the  pitching  moment  is  negative  (stabilizing)  in  constrast 
to  the  positive  value  for  the  ONERA  wing. 

The  Mach  number  contour  plot  of  Fig.  20  shows  a  shock 
pattern  on  the  surface  for  this  wing,  which  is  swept  forward 
near  the  tip.  and  this  shock  wave  is  stronger  than  that  ob¬ 
served  for  the  ONERA  wing.  Following  the  strong  shock 
wave,  there  is  a  massive  separation  of  the  flow  and  this  flow 
reattaches  only  at  about  the  trailing  edge  of  the  wing.  Also, 
the  shock  wave  extends  to  the  root  of  the  wing  as  before. 
The  shock  gets  weaker  towards  the  tip  region  and  the  flow 
stays  attached  in  this  region.  Figure  21  shows  the  surface 
oil  flow  pattern  for  this  wing.  The  lines  of  separation,  reat- 
tachment  and  other  important  flow  features  are  delineated 
in  this  figure.  Figure  22  shows  the  far-fteld  view  of  the  tip 
vortex  for  this  wing.  Also  shown  in  this  are  the  vorticity 
contours  in  the  Y-Z  plane  at  several  X-locations  of  the  tip 
vortex  behind  the  wing.  In  contrast  to  the  swept-tip  ON¬ 
ERA  wing,  the  tip  vortex  for  this  straight-  and  tapered-tip 
wing  does  not  lift  off  the  surface  of  the  wing  until  after  it 
passes  the  trailing  edge  of  the  wing  tip. 

The  nondimensional  tip  vortex  strength  determined  from 
integrating  the  velocity  vector  around  a  closed  path  sur¬ 
rounding  the  vortex  y  ields  a  value  of  0.17  w  hen  the  path  of 
integral  includes  the  vortex  sheet  of  the  wake  and  a  value  of 
0.21  for  the  path  excluding  most  of  the  wake  vortex  sheet. 
This  suggests  a  tip  vortex  strength  of  about  0.21. 

Conclusions 

A  multi  block  zonal  algorithm  for  solving  the  Euler  and 
the  Navier-Stokes  equations  is  used  to  numerically  simulate 
the  formation  and  roll-up  of  the  tip  vortices  of  wings  in 
subsonic  and  transonic  flows.  In  all.  four  different  wing 
geometries  have  been  used  as  test  cases.  The  influence  of 
the  tip-planform.  the  tip-cap  shape,  and  the  free-streatn 
Mach  number  on  the  formation  process  has  been  studied. 
The  numerical  result-  from  the  TNS  code  agreed  very  well 
with  the  results  from  a  single-bhx  k  PEN3I)  code.  Typical 
CPF  run-time  for  each  steady-state  solution  was  about  3 
hours  on  the  <  R  VY-XMP  machine. 

Comparison  of  the  numerical  result-  with  the  limited 
available  experimental  data  for  both  the  subsonic  and  the 
transonic  conditions  shows  good  agreement  for  the  surface 


pressures,  except  in  the  immediate  vicinity  of  the  tip,  and 
for  the  tip  vortex  strength,  where  available.  The  disagree¬ 
ment  of  surface  pressure  distributions  in  the  tip  region  is 
attributed  to  inadequate  capability  of  predicting  the  shock- 
induced  boundary-layer  separation  in  the  transonic  cases. 
The  subsonic  calculations  reproduce  the  qualitative  behav¬ 
ior  of  the  experimental  tip  vortex  formation,  including  the 
changes  in  the  tip  separation  and  in  the  vortex  lift  off  that 
have  been  observed  experimentally  and  which  have  been 
caused  by  rounded  tip-caps.  However,  the  square-tip  simu¬ 
lations  fail  to  produce  the  correct  suction  peaks  under  the 
vortex.  This  may  be  due  to  an  unrealistic  sharp  corner 
that  was  created  at  the  tip  leading  edge  by  the  H-H  grid 
topology.  Nevertheless,  it  has  been  demonstrated  that  it  is 
possible  to  calculate  the  initial  development  of  the  tip  vor¬ 
tex.  The  limited  study  presented  here  on  the  tip-cap  effect 
on  the  flow  field  near  the  tip  region  suggests  that  the  tip 
shape  is  an  important  ingredient  of  the  problem. 

As  expected,  the  formation  of  the  tip  vortex  involves 
braiding  of  the  particle  paths  in  the  tip  region  from  both 
the  upper  and  lower  surfaces  of  the  wing.  For  a  lifting 
wing  the  particles  from  underneath  the  wing  (high-pressure 
side)  cross  over  the  tip  region  on  to  the  upper  surface  (low- 
pressure  side)  and  modify  the  flow  field  in  the  tip  region 
by  the  three-dimensional  tip-relief  action.  The  tip  vortex 
first  lifts  off  from  the  surface  and  then  rolls-up  and  moves 
inboard  of  the  tip,  staying  distinctly  above  the  wake  vor¬ 
tex  sheet.  Both  tip-cap  modification  and  the  tip-planform 
change  had  similar  influences  in  increasing  the  gross  lift  of 
the  wing  and  hence  of  the  tip  vortex  strength.  In  addition, 
the  modification  to  the  ONERA  planform  delayed  the  lift 
off  of  the  tip  vortex  from  the  surface. 

The  results  demonstrate  that  although  there  is  still  room 
for  improvement,  realistic  three-dimensional  calculations  of 
viscous  flows  over  w  ings  and  their  associated  tip  vortex  for¬ 
mation  are  now  feasible.  This  computational  fluid  dynam¬ 
ics  capability  provides  a  new  tool  for  analyzing  and  im¬ 
proving  the  aerodynamic  characteristics  of  wings  and  other 
lifting  surfaces. 
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Fig.  3.  Surface  pressure  distributions  for  several  span- 
wise  stations  and  comparison  with  experimental  data  (Refs. 
16.17)  for  NACA  0015  wing.  -  0.16.  a  —  11  degrees, 
and  Re  =  2  million. 


Fig.  4.  Spanwise  cross  sectional  views  for  NACA  0015  wing 
showing  the  details  of  tip-cap  shapes  at  X  -  .5,  .75  and  .9. 
The  associated  grid  for  X  =  .5  station  shows  the  grid-point 
distribution  for  the  tip-cap. 


Fig.  6.  Surface  oil  fl 
0015  wing.  Mqo  —  0- 
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Fig.  7.  Close-up  view  of  the  surface  oil  pattern  in  the  tip 
region  and  the  lift  off  of  the  tip  vortex  for  the  NACA  0015 
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Fig.  9.  Surface  pressure  distributions  for  several  spanwise 
stations  and  comparison  with  experimental  data  (Ref.  21) 
for  Wing  C.  Moo  =  0.90,  a  =  5  degrees,  and  Re  =  6.8 
million. 
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Fig.  13.  Close-up  view  of  the  tip  vortex  for  Wing  C, 


Fig.  M.  Farfield  view  of  the  tip  vortex  for  Wing  C.  The 
vortirity  rountours  show  the  magnitudes  and  the  shapes  of 
the  tip  vortex  at  several  X-stations 
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Fig.  15.  Surface  pressure  distributions  for  the  ONERA 
wing.  V/^  0.85.  o  5  degrees,  and  Re  8.5  million. 
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Fig.  16.  Planform  view  of  the  upper  surface  Mach  number  Fig.  17.  Surface  oil  flow  pattern  for  ONERA  wing  showing 
contours  for  the  ONERA  wing.  separation  and  reattachment  lines. 
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Fig.  20.  Planform  view  of  the  upper  surface  Mach  number  Fig.  21.  Surface  oil  flow  pattern  for  the  modified  ONERA 

contours  for  the  modified  ONERA  wing.  =  0.85,  a  =  wing  showing  the  separation  and  reattachment  lines. 

5  degrees,  and  Re  =  8.5  million. 


